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TIME -DEPENDENT STRUCTURE OF ONE -DIMENSIONAL MAGNETIC 


COMPRESSION WAVES IN COLLISIONLESS PLASMAS - 
OBLIQUE AMBIENT MAGNETIC FIELD 
By Vernon J. Rossow 
Ames Research Center 

SUMMARY 


Numerical analysis is made of the structure of one -dimensional unsteady 
magnetic compression waves propagating into collisionless plasmas when the 
ambient magnetic field is oblique to the wave direction. In the cases studied, 
both the subcritical and supercritical solutions changed continuously from one 
limiting situation to the other. A similarity rule was found for correlating 
the subcritical cases so that the two -parameter family of solutions in mass 
ratio and field angle were reduced to a single set. 


INTRODUCTION 


Recent studies of compression waves propagating through collisionless, 
completely ionized gases indicate that a large variety of wave shapes is 
possible. Such analyses bring about an understanding of these waves, and 
hopefully, a compilation of the many possible situations will help the experi- 
menter identify the origin of probe data obtained in space or in the labora- 
tory. With such considerations as motivation, the present study was carried 
out to find the structure of one -dimensional unsteady magnetic compression 
waves moving into a two -component, completely ionized plasma that is, in turn, 
embedded in a magnetic field. The initially cold plasma is compressed by a 
magnetic field that is driven by an electric field (the left boundary in 
fig. l) . The magnetic field in the undisturbed plasma has two constant compo- 
nents B xo and B zo , with tan 0 = B zo /B xo . It is assumed that the flow field 
can be rotated about the x axis to eliminate any y component of the ambi- 
ent magnetic field. The disturbance field, Ey, at the initiation plane causes 
the magnetic piston to be composed of only B z and thereby imposes a restric- 
tion on the possible orientations between ambient and disturbance fields. In 
general, the compressed layer will involve all three components of the mag- 
netic field and of the velocity vector. Of particular interest here is the 
variation in the structure of the compressed layer as the angle 0 between 
the direction of wave propagation and the ambient magnetic field is changed 
from the alined (0 = 0°) to the transverse (0 = 90°) situation. 

The structure of a single or solitary wave has been treated by Adlam and 
Allen (ref. 1) for a wave traveling perpendicularly to the ambient magnetic 
field and by Saffman (ref. 2 ) and Pataraya (ref. 3) for a wave moving along 


the field lines. The results for the two limiting cases were generalized by 
Saffman (ref. 4 ) and Kellogg (ref. 5 ) to the cases where the propagation direc- 
tion is at an arbitrary angle to the ambient magnetic field. Rather than just 
a single pulse , the general case consists of a group of pulses that may or may 
not be stable (refs. 4 and 5 )- Morton (ref. 6) studied the development with 
time of a one -dimensional compressed layer brought about by a magnetic piston 
moving at any angle relative to the undisturbed magnetic field direction. He 
found by an approximate treatment that the structure of the wave changes 
rather abruptly near the transverse limit when the mass ratio R = mi/me is 
large. For stronger waves, his more exact results obtained by numerical inte- 
gration of the equations indicate that the change was not so abrupt as first 
thought; consequently, results near 9 = 90° resemble those at 90° • 

The present investigation first supplements the work of Morton (ref. 6) 
by redoing the so-called subcritical waves 1 by the method of Auer, Hurwitz, 
and Kilb (refs. 7 and 8). Morton’s results are then extended by new results, 
found by the present technique, for the supercritical case in which breaking 
occurs and wherein volume elements of plasma interchange positions during the 
event on a scale comparable with the size of the compressed layer. Both the 
present analysis and that of Morton start with the same basic equations but 
differ in the manner of solution. Whereas Morton used a continuum formulation 
with an artificial viscosity in the numerical work, here Lagrangian coordi- 
nates are used to follow the plasma elements in order to integrate the time- 
dependent flow field. Such a discrete model permits the calculation to 
proceed through breaking of the waves and interchanging of particles without 
requiring more assumptions as to the nature of the process. 


DEVELOPMENT OF EQUATIONS FOR NUMERICAL ANALYSIS 


Introducing the longitudinal (B xo ) in addition to the transverse 
component (B zo ) of the magnetic field causes the three velocity components to 
be coupled to each other through the electromagnetic fields. As a result, the 
procedure used here differs from that of Auer, Hurwitz, and Kilb (refs. 7 
and 8) only in that more components of the various quantities must be carried 
along in the analysis. A brief description of the equations leading to the 
numerical work will be given in this section and emphasis will be given to the 
new material. 

In vector notation, the equations for conservation of mass and momentum 
of the ions (subscript i) and of the electrons (subscript e) for the planar 
flow field under consideration are 

■‘■Convention has defined subcritical compression waves as those in which 
the plasma elements retain their original order. By contrast, supercritical 
waves are those wherein the wave intensity is so large that plasma elements 
exchange positions on a large scale so that the structure of the compressed 
layer is irregular. When 9 = 90 ° , the division between the two types of 
waves occurs at an Alfv£n Mach number of 2.0; 

l/2 

Alfven Mach number = M A = U wave /{B 0 /[pn 0 (m e + mi)] 7 } 
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since all the derivatives with respect to y and z are zero and the ions are 
singly ionized. The quantities nq^ e , mq^ e are the number density and mass 
of the ions and electrons, respectively, and the velocity is U = iu + jv + few. 
Similarly, the electromagnetic field equations for the one -dimensional wave 
system are 
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where the displacement current has been assumed to be negligible. Since the 
plasma frequency is many orders of magnitude greater than any of the cyclotron 
frequencies, the number density of ions and electrons are taken to be equal as 
a first approximation (ref. 9) • Therefore, 

nq = n e = n and uq = u e = u ( 7 ) 

As in previous analyses (refs. 7 and 8), the electrons are assumed to retain 
their original order while the ions are permitted to take any order that their 
dynamics might dictate. Justification for this approximation is given in 
reference 9. 


As an aid in the numerical work, the magnetic vector potential, A, and 
the electric potential, <p, are introduced such that 

dA 

B = V X A and E = - ^7 - VCp 

For the one -dimensional unsteady problem being treated here. 
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These equations can be integrated to yield 


Ay(x,t) = 

f B z (x 1 ,t)dx 1 - r 
0 '■'o 

' Ey(0,t 1 )dt 1 + Cy 

(9a) 

A z (x,t) = 

( x lj "t) ^1 

J*° E z (0,t 1 )dt 1 + c z + yB xo 

(9b) 


where the constants C y and C z are determined by conditions at the reflection 
or right side of the flow field. The quantity yB xo is added to complete the 
definition of the magnetic vector potential and does not enter the analysis 
because the y coordinate can be ignored. 


Before the differential equations can be further adapted to their 
difference form for the numerical calculations, the initial conditions (i.e. , 
at t = 0) and the restrictions at the boundaries x = 0 and x = x m need to 
be specified. As assumed by Auer, Hurwitz, and Kilb (refs. 7 and 8) and By 
Morton (ref. 6), the particles are taken to be cold and therefore stationary 
at t = 0 with uniform distribution throughout the flow field. The magnetic 
field in the undisturbed plasma has the two constant components B xo and B z0 
(see fig. l) . At the initiation plane, the disturbance to the plasma is 
assumed to be brought about by an electric field in the y direction 

E y (0,t) = E yo (l - e" at ) ; E z (0,t) = 0 (lO) 

that generates a magnetic field in the z direction. (The time-rise factor 
a was always taken to be 0.12 in the calculation.) This magnetic field 
forces the plasma away from the left boundary with a piston-type action so 
that a wave field and compressed layer are produced. It is assumed that the 
reflection plane at x = x m is impermeable to both the magnetic field and the 
plasma so that all disturbances are reflected without loss. The magnetic vec- 
tor potentials there are taken as Ay(x m ,t) = A z (x m ,t) =0 so that 
Cy = -x rft B z0 and C z = x m By 0 = 0 in equations (9) • 


With these quantities as starting conditions, the y and z momentum 
equations for the ions and electrons can be integrated along the particle 
paths to yield 


fit . 


eA. 


o+ 

y 


eB. 


’xo 


W.? 


w e ' 


mi 

eA. 

eA z 

“T 

eA z 

me 


mi 


m. 


eA. 


o- 


i 

eB. 


xo 




“e 


eA 


o+ 


eB. 


xo 


mi 


mi 


eA, 


rn 


o - 


eB. 


xo 


me 


(Ha) 

(lib) 

(He) 

(lid) 


4 



where the velocities are zero at t = 0 and the superscripts on A designate 
the value of the magnetic vector potential at the location of the particle at 
t = 0. 
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Capital letters are now used to designate the variables in dimensionless 
form where the reference length, time, velocity, etc., are given by 


* = U A/ U ch 

t* = a/u* 

U* = (u A u D ) l/2 

V* = W* = (mj_ + m e )U*/(mj_m e ) 172 
B* = B q U*/U a 
A* = AB* 

with B 0 = (Bxo + Bf 0 ) l/2 , U a = Bo/tprioCmi + me)] 172 , U D = E y0 /B 0 , and 
u> ck = eB 0 /(mj_me) l72 . The dimensionless variables then become U = u/u*, 

BY = By/B*, AY = Ay/A*, X = x/A, T = t/t*, etc., and the mass ratio is defined 
as R = mi/me . It is also convenient to introduce the Lagrangian variable, £, 
so that 


I = 




As was done previously, the flow field is divided into slabs of plasma, 
the thicknesses of which are then made vanishingly small. Each slab or sheet 
contains the same number of ions and electrons and is located at the mass 
center of the element it represents. The procedure used in reference 8 is 
used to set up an array of equations to find the motion of the slabs of plasma 
and the associated electromagnetic fields. Recursion relations used in the 
machine program to find the magnetic field components throughout the flow 
field are then found for the kth slab as (compare, e.g., with ref. 10) 

BY k _ 1 - [2 + D(X k+1 - X k )]BY k + BY k+1 + HY k = 0 (l2a) 

BZ k-i ' t 2 + D ( x k+i " x k)l BZ k + BZ k+i + HZ k = 0 (12b) 

where D = X^N and 

HT k = ^k+i ' XB k + R ( YE k+i ~ ra k )] (13a) 

HZt - ~ ~ U b^ 1 ~ ~ + [ZI k +1 - ZI k + R ( ZB k + i - ZE k)] 

(13b) 
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since BYO = 0. This matrix is inverted for the magnetic field values by 
first calculating the quantities 


RY k = l/[2 + D(X k+1 - X k ) - RY k+1 ] 

(l4a) * 

SY k = RY k (SY k+1 + HY k ) 

(Ub) ^ 

RZj, = l/[2 + D(X k+1 - X k ) - RZk +1 ] 

(14c) 

SZ k = RZ k (SZ k+1 + H^) 

(I4d) 


starting at the reflection plane where k = N and proceeding to the initia- 
tion plane magnetic piston region where k = 1. Since the boundary at X m is 
a perfect reflector, the beginning values for k = N are 

RY n = 1/[1 + 2D(X m - X N )] 

sy n = 0 

RZfj = 1/[1 + 2D(X m - Xjj) ] 

SZjj = RZjj ( HZjj ) 

Once these quantities are known, the recursion equations 

BY k = RY k (BY k _J + SY k (l5a) 

BZj^ = RZ^BZ^) + SZ k (15b) 

yield the local magnetic field after the value at the initiation plane (and in 
the piston region) is determined by 


BY™ = 


SYi - (BX0 /R i/2 )[r(YE 1 ) + YIt ] 
1 - RY X + D(Xi) 


(l6a) 


SZn + D i -AY,, + 


BZw - 


P(BZO) 


w R + 1 
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- N(R + 1) + 


AY 


o+' 


BZO 


+ [r(zei) + zin 


1 - RZi + D(Xi) 


(l6b) 

where AY W = AY(0,T) = -[T + (e^ - l) /a] - X m (B zo ) , since EYO = 1. The 
acceleration of the slabs of plasma is found by calculating the magnetic 
pressure difference on its two sides; that is, 

x k = -(BY* - BYg_ x + BZ^ - BZ|_ 1 )/2D (lT) 

The accuracy of the numerical calculations can be monitored to a certain 
extent by checking to see whether the momentum and energy content of the 
entire flow field is conserved. The three components of the momentum are 
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evaluated by means of 


2D?U k = r T [BYW 2 + BZW 2 - BY N 2 - BZ N 2 ]dT 3 
1 °o 


N , , /n N, 

2[VE k + R(vi k )] = R 1/2 BXO S(ZI k - ZE k ) 

|[WE k + R(WI k )] = -R l/2 BXO |(YI k - YE k ) 


(l8a) 

(l8b) 

(18c) 


An energy check during the calculations is made by evaluating and comparing 
the two sides of the expression 


iji N 

2 J (BZ w )EY(0,T 1 )dT 1 = D^ U k 2 + (R + 1) (vi k 2 + WI k 2 + 


VE k 2 + WE k 2 


+ ( X k+l ' X k)( BY k 2 + BZ k 2 ) 


( 19 ) 


A machine program for solving the foregoing system of equations was used 
to find a number of solutions. Results from the computer were written on 
magnetic tape and then plotted electrically. Accuracy tests were carried out 
on various cases to determine which mesh sizes were most suitable. 


RESULTS OF NUMERICAL ANALYSIS 


Emphasis is first placed on the variation in the structure of the 
compressed layer as the angle 0 between the wave direction and the ambient 
magnetic field changes from alined (0°) to normal (90°) incidence. The effect 
of mass ratio on the solution is then considered for two situations: first, 

the ambient magnetic field is chosen to be 1.0 so that the compression wave is 
subcritical and the plasma sheets retain their original order, that is, 

B q = (BXO 2 + BZ0 2 ) 1/2 = l.o(M A 1-5); second, the structure is considered for 
B 0 = 0.2 (M a 5.8). For a wave of this strength, the plasma slabs interchange 

order on a large scale and a random irregular -type flow field results. All 
other factors in the calculations are held fixed as 0 varies from 0° to 90 . 


Subcritical Case 

Data on the structure of a number of the subcritical cases are presented 
in figures 2 through 10. Of particular interest here is the manner in which 
the structure changes as the magnetic field angle, 0, or mass ratio, R, is 
varied. The growth of the compressed layer with time is shown in figures h 
through 10 for selected cases so that the reader can see how the flow field 
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shown in figures 2 and 3 developed. Variation of the velocity components with 
distance at a given time are presented in part (a) of each figure along with 
the magnetic field components. Special features of each solution are noted in 
the figure titles. % 

Figure 2 presents the structure of the compressed layer for & 1.5 at 
several angles when the time T is 60, that is, 60 time units after the dis- I 
turbance was initiated. At the mass ratio chosen, the flow field changes 
gradually from one type to another - the changeover appears to occur at about 
60°. Although the calculations were not made specifically for comparison, the 
curves for the magnetic field intensities are in general agreement with the 
solutions of Kellogg (ref. 5) and Morton (ref. 6) and no striking disparities 
are noted. 

Several changes in the flow field are observed as the angle of the 
ambient field varies. Because of the presence of a precursor -type wave, the 
compressed layer or disturbed region is thickest at 0=0° (see also 
figs. 5-9) • Also apparent in figure 2 is the extent that the magnetic field 
associated with the disturbance diffuses into the compressed layer at 0=0° 
causing the piston region to be smaller there. In contrast, this mixing at 
the piston-compressed layer interface is only several units of A at 90 ° • 

Although most of the structures shown in figures k through 10 appear to 
have established a pattern, they may not have been fully developed when the 
calculations were terminated. (Machine storage capacity prevented an increase 
in the size of the flow field and, hence, in the running time. Had the pres- 
ent cases been carried further, excessive interference from the reflection 
plane would have resulted.) The rather abrupt appearance of the small wave 
sequence in the 0 = 60° case (see fig. 8(b)), at about T = 35> suggests that 
prolonged tests might have brought a new wave system to light in some of the 
other solutions. This does not seem likely though because the structure of 
the compressed layers appears to have settled down to a pattern. 


Similarity Parameter for Subcritical Flow Fields 

Whereas the subcritical cases for 0 = 90° are independent of the mass 
ratio R, all other angles are affected by it. Both Morton and Kellogg esti- 
mated that the flow field changes rapidly with angle near 90 ° when the mass 
ratio is large. In order to study the effect of R on the structure of the 
compressed layer, a sequence of runs was made for 0 = 60 ° and the results at 
T = 60 are presented in figure 3* Although nothing new appeared the structure 
changed considerably. Of most importance, however, is the resemblance between 
the three new solutions and the results for 0 = 90 ° > 75° > and ^4-5° in figure 2 
(the 60 ° cases for R = 25 are duplicates). These results led the author to 
find an empirical rule for relating one case to another and then to the justi- 
fication of the results by means of the equations. A best fit to the data in 
figures 2 and 3 was obtained with the rule 

BZOj = BZOjj 
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bxo i 2 (r i 1/2 - 1) = bxo ii 2 (r ii 


1/2 


- 1 ) 


where I and II are the two cases being related by the rule. However, an 
examination of the recursion equations and, in particular, the quantities HI 
, and HZ (eqs. (13)) led to the similarity parameter that is felt to be the 
correct one. 


BZOj — BZOjj 


(20a) 


BXO-j-^R-j- 172 - Rj' 1 / 2 ) = BXOj^CRjj 172 - Rjx" 172 ) 


(20b) 


The two forms are equal at R = 1 and become equivalent at large mass ratios 
since R l72 - 1 = (R - l)/(R 172 + l) « (R - l)/R 172 . Equation (20a) arises 
from the fact that the solutions at 90° are independent of R as long as the 
plasma slabs do not interchange order. Equation (20b) is justified by consid- 
ering the lateral displacements of the particles in the two cases as they 
appear in equations (13). Since the mass of the electron is the same for both 
cases, the relative lateral displacement is reduced by an increase in BXO, so 
that, for example. 


Z e SK -const ant/BXO 

Lateral movement of the ion is approximated by 

Z^ ex +(R constant) (l/R) = constant/BXO 

The first term describes the effect on the Larmor radius brought about by the 
longitudinal magnetic field and mass of the ion and the second term treats the 
reduced acceleration (and, hence, lower velocity and distance) for a heavier 
ion. Hence, the ion and electron have the same proportionality factor but 
have opposite signs. When these results are inserted into the last term of 
equations (13), which involves BXO/R 172 , the similarity rule presented in 
equations (20) is derived. 

As a test of the similarity rule (20), two sets of related cases were 
calculated and the results are compared in figure 4 at three time intervals. 
Exact correspondence is not achieved, but the general features of the magnetic 
field lines are reproduced quite well considering the large difference in mass 
ratio between the cases compared. As is to be expected, the other flow-field 
parameters deviate more from each other. A horizontal line for the piston 
region on BZ was omitted to illustrate the difference in the motion of the 
interface. It appears that the average velocity is about the same for the 
cases compared but the heavier ions oscillate at a lower frequency and greater 
amplitude than the lighter ions. The high frequency oscillations of the 
lighter ion are not duplicated (fig. 5(b)) by the heavier ion. As a result, 
the magnetic field of case III (R = 1836 ) appears to be a mean of case IV 
(R = 25). From these results, the similarity mile, equations (20) , appears to 
give satisfactory similar solutions. 
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Supercritical Flow Fields 


The variation of the magnetic field components in the compressed layer 
for a disturbance of supercritical strength is first illustrated in figure 11 
for several values of 0 from 0° to 90°. Once again, the structure does not 
change abruptly from one form to the other. Transition from the 0° pattern to 
the type typical of the transverse situation occurs at about 30°. The forma- 
tion of a region of nearly uniform magnetic field near the piston also appears 
to be characteristic of the solutions near 0°. As 0 increases from about 
30°, the extent of this area decreases quickly to zero. Since ®^)piston an ^ 

the length of the piston region are essentially independent of 0, the uniform 
compressed layer near the piston for 0=0° and 15 ° does not appear to be 
caused by diffusion of the disturbance magnetic field into the compressed 
layer. It may be the beginning of the formation of a uniform compressed layer, 
in which case, the forward part of the layer could be considered as the shock 
structure (in a gasdynamic continuum sense). The calculations could not be 
carried further to see if this were true. Whether the depth of the compressed 
layer is a function of 0 cannot be judged from figure 11 because of the 
unsteady nature of the supercritical cases. Disturbances at X = 140 to l60 
for the 0 = 15 ° solution result from particles thrown far ahead of the com- 
pressed layer by the start-up of the piston field. Forward spraying of the 
particles occurs in most of the strong or supercritical solutions but such 
plasma motion is not always as apparent in the magnetic field profiles as 
exhibited for 0 = 15° in figure 11. Although the z component of the mag- 
netic field exhibits a random profile of about the same magnitude at all 
angles, the intensity of the y component decreases regularly from 0° to 90 °. 
For angles of 4-5° and above, the magnitude of BY relative to BZ appears to 
be the most obvious indication of the obliqueness of the ambient field. 

The influence of the mass ratio on the wave structure is shown in 
figure 12 for 0 = 60°. This angle was chosen to parallel the subcritical 
solutions. Note that the magnetic profiles do not change with R in quite so 
pronounced a fashion as did the subcritical cases in figure 3. A slight 
increase in the size of the disturbed region for R = 100 is attributed to the 
unsteady character of the flow field and not to an increased average growth 
rate. These solutions have the same irregular shape as those in figure 11, 
but they do not resemble any of them closely enough to form a correlation. In 
particular, the time and distance dilation formula derived empirically by 
Auer, Hurwitz, and Kilb (ref. 8) for the 90° supercritical waves does not 
apply for any other angle because the depth of the compressed layer does not 
change with R. Attempts were made to arrive at some sort of similarity rule 
for these results, but no satisfactory parameter was found. 

As was done for the subcritical cases, the development with time is 
shown for several cases in figures 13 through 17 • Special features of the 
solutions are again noted in the figure titles. 
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GONG HIDING REMARKS 


Calculations of the flow field generated in a collisionless plasma by a 
disturbance of the magnetic piston type yielded results that agree in general 
with previous work on oblique magnetic compression waves. The structure 
changed continuously with angle from one limiting form to the other, that is, 
from 0° to 90°. Although the solutions at intermediate values of 0 are more 
complex in their dependence on mass ratio than the results at 90°, a similar- 
ity rule was found for subcritical wave speeds. This rule makes it possible 
to relate magnetic field results at one mass ratio and angle to a correspond- 
ing one, thereby reducing a two -parameter family to a single set. A compara- 
ble rule was not found for the supercritical flow fields. The need for such 
a rule is not so great because the solutions do not depend strongly on mass 
ratio . 

Limitations imposed by computer capacity did not permit the calculations 
to be carried out as far as desirable in all cases. Although it does not seem 
likely, some of the foregoing conclusions might have been different had the 
flow field been allowed to develop for a sufficiently longer period. Also, it 
should be remembered that the disturbance used to generate the wave was an 
electric field in the y direction at the initiation plane. Such a field 
produces a magnetic piston composed of B z only. Another type of piston 
might yield quite different results. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., 94035, Feb. 2, 1967 
129-02-03-03-00-21 } 
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Comparison of two test cases as predicted by similarity rule. 
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(a) Flow-field variables at T = 60. (Scatter in the data at the piston- 
plasma interface is characteristic of solutions near 0=0° and appears 
to result from an accumulation of plasma slabs at that location.) 

Figure 5.- Structure of time -dependent flow field for subcritical wave; 

R = 25, 0 = 0°, B c = 1.0 (BX0 = 1.0, BZO = 0) . 
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(a) Flow-field variables at T = 60 . (About the only difference between this 
case and the one at 0 = 0° is the increased area under the BZ curve 
brought about by its present nonzero value in the ambient field.) 

Figure 6.- Structure of time -dependent flow field for subcritical wave; 

R = 25 , 0 = 30°, B q = 1.0 (BXO = 0 . 866 , BZO = 0-500). 
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Figure 7 .- Structure of time -dependent flow field for subcritical wave; R = 25, 0 = 45 °, B Q = 1.0 
(BXO = O. 707 , BZO = 0 . 707 )- (Profiles of the flow variables at T = 60 are not available because 
of difficulties with the computer program when this case was calculated. These results are 
presented because the BY field appears to be developing a sequence of separate (or solitary) 
wave trains similar to those discussed by Kellogg (ref. 5)-) 




(b) Development of magnetic field components with time. (As with the 9 = 45° case, separate 
trains are developing in BY. Now, however, the more fully developed set that is out for 
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Figure 8.- Concluded. 
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(a) Flow-field variables at T = 60. (Results for this case closely resemble 

those for 0 -- 90° (ref. 10).) 


Figure 9- - Structure of time -dependent flow field for subcritical wave; 
R = 25, 0 = 85°, B 0 = 1.0 (BXO = O.O872, BZO = 0.9962) . 






Figure 10 .- Structure of time -dependent flow field for subcritical wave; 
(BXO = 0 . 500 , BZO = 0.866). (Note that BY = 0 , BZ is the same as a 
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(a) Flow-field variables at T = 60. (Points near U ** 2 in the upper curve 

are plasma sheets sprayed forward by the rapid initiation of the distur- 
bance field. Such motion of the plasma in the X direction is not 
inhibited because the transverse magnetic field is zero,) 

Figure 13 .- Structure of time -dependent flow field for supercritical wave; 

R = 25, 0 = 0°, B q = 0.2 (BXO = 0.2, BZO = 0) . 
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(a) Flow-field variables at T = 120. (A thermalized layer near the piston 
also builds up in this case. Forward spraying of slabs by the distur- 
bance field at T ^ 0 is reduced by the nonzero value of BZO. Waves in 
the compressed layer appear to not resemble the solitary waves of Kellogg 
(ref. 5) so much as the waves typical of 0 = 90 ° (ref. 10) . ) 

Figure 14-.- Structure of time -dependent flow field for supercritical wave; 

R = 25, 0 = 30°, B q = 0.2 (BX0 = 0.1732, BZO = 0.100). 
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(a) Flow-field variables at T = 120. (Although the compressed layer is 
nicely thermalized, this solution appears not to reduce to that at 
0 = 90 ° as in the subcritical case.) 


Figure 1 6 .- Structure of time -dependent flow field for supercritical wave 
R = 1, 0 = 60°, B 0 = 0.2 (BXO = 0.100, BZ0 = 0.17^2). 









*'The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 
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